هی ی 9 101:1022067/67332 سال سی و منوم. شماره دیء ۱۳۹۹ 


تحلیل بایداری یک تونل در شرابط هیدروستاتیک با استفاده از روش‌های مختلف قابلیت اعتماد در 
نرم‌افزارهای 1۲ و »اونع۵1* 
مقاله علمی - پژوهشی 


هادی فتاحی(٩‏ فاطمه جیریایی شراهی( 


چکیده یکی از راه‌های تحلیل پایداری تونل, بررسی همگرایی و شعاع منطقةٌ پلاستیک در اطراف آن است. از طرفی وجود عدم قطعیت در 
باراظرهای طراصی :ای ایعتی شبازه‌ها ره اسفاهه نرق سای امالاتی و قابلیت اضتماه متری مر دف رف آین تین متعی الیه ی 
بهره‌گیری از روابط هوک و براون و موه رکولمب. یک متدولوژی برای تحلیل پایداری دو توئل (در شرایط متفاوت) با | ستفاده از روش‌های 
قدرتمند قابلیت اعتماد شامل روش مرتبهُ اوّل (/0) روش مرب دوم )٩010(‏ و شبیه‌سازی مون تکارلو در نرم‌افزارهای جدید 111و 
68 ارائه شود. پس از مدل سازی احتمالاتی در اين دو نرم‌افزار, نتایج تحلیل‌ها ذ شان می‌دهد که روش قابلیت اعتماد مرت دوم ذ سبت به 
روش مرتبهُ اول در توابع پیچیده‌تر حالت حدی, از دقّت بیشتری برخحوردار است. همچنین در مقادیر بالای احتمال شکست, روش‌های مرتبه 
اول و شبیه‌سازی موذ تکارلو حروجی‌های واقعی‌تری را ارائه می‌دهند و در مقادیر پایین احتمال شکست» روش شبیه‌سازی مون تکارلو توصیه 
نمی شود. به‌علاوه تونل مدنظر در شرایط بدون نگهداری دارای احتمال شکست بالایی (2۲6) است و تدارک یک سیستم نگهداری با فشار 
نگهداری ۰/۱۲ مگاپاسکال کارامد است. 


واژه‌های کلیدی قابلیت اعتماد نرم‌افزار 11 نرم‌افزار 61856 تحلیل پایداری تونل. 


هام )۱۱۲/۵۵ منوا عصمتانلصم اهاو00 ۲۳۱۲ جصز تمصصت 1 و ۵۲ عتعراممم وانانماهاد 
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۰ 5۲2011107 [مصصها [ رنه اگمو اون ۵ رمیهسامو ]۳ روانازماهنای۲ ۷۷۲۵۲۵ مک 


* تاریخ دریافت مقاله ۹۸/۱۲/۲۰ و تاریخ پذیرش آن ۹۹/۵/۱۳ می‌باشد. 
(۱) نویسندة مسئول دانشیار دانشکدة مهندسی علوم زمین. دانشگاه صنعتی اراک. رم انقصصع ۵ 62 تطم)نع! :انهصصرز 


(۲) کارشناسی ارشد. دانشکده مهندسی علوم زمین دانشگاه صنعتی اراک. 


مقدّمه 
تونل با توجه به وسعت کاربردش در کارهای عمرانی و 
معدنی از دیرباز نقش مهمی در تسهیل فعالیت‌های بشر 
داشته است. لذا در سال‌های اخیر احداث تونل‌ها به‌طور 
چشمگیری افزایش یافته است [1-3]. حفر فضاهای 
حفرشده می‌شود و منطقةٌ پلاستیکی با پتانسیل تخریب 
در اطراف این فضا ایجاد می‌کند [5 ,4]. شرط ادامهٌ کار 
يا استفاده از این فضاء وابسته به پایداری و ثبات آن در 
و بررسی شعاع منطقهٌ پلاستیک و پیش‌بینی عامل 
نگه‌دارنده برای آن امری لازم است. از طرفی همگن یا 
همسان گرد بودن خاک و سنگ. یک خاصیت ایده‌آل برای 
آن است که هیچ گاه تحقق نمی‌یابد و ویژگی‌های آن بسته 
به ضریب تغییراتش از نقطه‌ای به نقطةٌ دیگر تغییر می‌کند. 
موجب پی بردن به رابطه یا همبستگی بین این ویژگی‌ها 
نیز می‌شود. در اصطلاح به این تغییرات در ویژگی‌ها.؛ 
عدم قطعیت گفته می‌شود. مهندسی. تصمیم گیری تحت 
شرایط عدم قطعیت است و انکار این نکته امکان‌ناپذیر 
است. قطعیت در متغیرهای مدل به معنی در نظر نگرفتن 
طراحی‌ها و مدل‌ها راهی برای ورود عدم قطعیت باز شود 
[5 ,4]. به‌منظور در نظر گرفتن عدم قطعیت در تحلیل‌های 
مهندسی احتمالاتی و قابلیت اعتماد 
که درنهایت با استفاده از روش‌های مختلف عددی و 


روش‌های 


شبیه‌سازی به دست می‌آید» همان احتمال شکست 
(وانلژطه۲:00 معسانه۳) يا احتمال خرابی است. احتمال 
شکست را می‌توان به شکل زير تعریف کرد [6]: 


لر -( > ۳60۵0 << 


2 (٩) 


در اين رابطه ۲۴ احتمال شکست. (8)16 تابع حالت حدی 


(ممناعص؟ عاماه انهن) سیستم مایت رین] یر 


نشريه مهندسی عمران فردوسی 


تحلیل پایداری یک تونل در شرایط هیدروستاتیک با استفاده از.. 


بردار متغیرهای تصادفی و (06] تابع چگالی احتمال 
شترک (۳۱۵08 اتوجهنا بانلامادامنج اصژم) همه 
متغیرهای تصادفی است. به عبارتی این انتگرال وظیفه 
دارد مساحت مشخص‌شده در شکل (۱) را به دست آورد 
[6]. 


#ین 


شکل (0۱: شاخحص قابلیت اعتماد 1 و احتمال شکست [7] 


حل این انتگرال در حالت عمومی بسیار پیچیده و 
غیر ممکن است؛ به اين دلیل که ممکن است تابع حالت 
حدی» یک ترکیب غیر خطی باشد و يا متغیرها توزیع 
نرمال نداشته باشند؛ بنابراین برای رسیدن به جواب. باید 
به راه‌های تقریبی روی آورد. ازاین‌رو با پذیرش چند 
فرض ساده‌کننده به‌جای محاسبة دقیق احتمال خرابی» از 
پاراتر ‏ دیگری ‏ ه نام شاخص قابلیت ‏ اعتماد 
ر :ععهع بوانانطاه‌نام0۳۴) استفاده می‌شود تا حجم 
فحاسات وودان مووو از کاهتی انش روش‌مای توق 
کب ایس قافن کرش کی ارت اند اد رو مرف 
اول (۳۵۲۴۷۲ :۱۲۵00 باتازمدتام؟ 0:۳0 افا) و 
مرتبهةً دوم قابلیت اعتماد ( 7انلطدنام؟ 0:06 ٩66020‏ 
7 :۷]60000)._ همچنین شبیه‌سازی مونت‌کارلو 
(۷۵ :ممتاماصنو ماه ۷]0868) روشی است که با 
داده‌سازی زیاده احتمال شکست را به دست می‌آورد. 

محققین بسیاری از دیرباز متوجه برخی عدم 
قطعیت‌ها در طراحی و ساخت سازه‌های ژئومکانیکی 
شده بودند. رویکرد برخی از آنها نسبت به عدم قطعیت؛ 
توسط روش‌های قابلیت اعتماد بدین شرح است: 


سال سی و سوم شماره یک, ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی 


سانگ ان چو در سال ۲۰۰۹ با ترکیب روش مرتبة 
دوم قابلیت اعتماد و شبکة عصبی مصنوعی به تحلیل 
احتمالاتی شیروانی خاکی پرداختند [7]. سانجی و کیچی 
در سال ۲۰۰۹ بر روی روانگرایی خاک بر اساس داده‌های 
آزمایش نفوذ استاندارد. آنالیز قابلیت اعتماد مربوطه را 
انجام دادند [8]. چین لونگ و بک در سال ۲۰۱۲ شمع‌ها 
یا المان‌های انتقال بارهای سطحی به عمق زمین را با 
استفاده از احتمالات و قابلیت اعتماد تحلیل کردند [9]. 
پنگ و همکاران در سال ۲۰۱۳ به تحلیل برخحی مسائل 
ژئوتکنیکی نظیر نشست فونداسیون سطحی. ظرفیت 
باربری یک شالودهُ سطحی و پایداری شیروانی خاکی با 
استفاده از روش‌های قابلیت اعتماد پرداختند [10]. 
لوائیس و فیلیپ در سال ۲۰۱۷ پایداری داخلی دیوار 
خاکی تقویت‌شده را با روش مونت‌کارلو مورد بررسی 
قرار دادند [11]. رخشانی‌مهر و همکاران در سال ۱۳۹۲ 
ایمنی شمع‌های بتنی و فولادی تحت بار ثقلی برای دو 
نوع خاک چسبنده و غیر چسبنده را با روش‌های 
شبیه‌سازی مونت‌کارلو شبیه‌سازی وزنی و مرتبة اوّل 
قابلیت اعتماد بررسی کردند [12]. هوک در سال ۱۹۹۸ 
آنالیز قابلیت اعتماد یک تونل دایره‌ای را به روش 
شبیه‌سازی مونت‌کارلو با استفاده از نرم‌افزار ریسک را 
انجام داد [13]. لی و لو در سال ۲۰۱۰ برای محاسبة 
شاخص قابلیت اعتماد یک تونل دایره‌ای تحت تنش 
هیدروستاتیک از روش مرتبهٌ اوّل قابلیت اعتماد استفاده 
کردند [14]. یانگ و همکاران آنالیز قابلیت اعتماد را برای 
بررسی پایداری سقف تونل‌های عمیق در دو لایه توده 
سنگ به کار بردند [15]. لو و اینستین در سال ۲۰۱۳ در 
مورد گوه‌های سقف و به‌کاربردن نیروی راک بولت در 
تونل آنالیز قابلیت اعتماد را برای اين تونل انجام دادند 
[16]. لی و لو [17]به کمک روش‌های قابلیت اعتماد 
مرتبةٌ اوّل و روش شبیه‌سازی مونت‌کارلو به تحلیل 
پایداری یک تونل با مقطع دایره‌ای تحت‌فشار 
هیدرواستاتیک پرداختند. لو و همکاران [18] به آنالیز 


سال سی و سوم شماره دی ۱۳۹۹ 
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قابلیت اعتماد اندرکنش بین سنگ و سیستم نگهداری در 
تونل‌های دایره‌ای با استفاده از روش‌های قابلیت اعتماد 
مرتبهٌ اوّل و سطح پاسخ پرداختند. لو و و [19] به کمک 
معیارهای هوک-براون و موهر-کولمب و با استفاده از 
روش سطح پاسخ و روش‌های قابلیت اعتماد مرتبة ال 
و دوم به بررسی شعاع منطقه پلاستیک و جابه‌جایی تونل 
پرداختند و درنهایت نتایج را با روش شبیه‌سازی 
مونت‌کارلو مقایسه کردند. گه و ژانگ [20] به کمک 
روش‌های قابلیت اعتماد ابتدا یک مدل رگرسیون 
لگاریتمی را برای پیش‌بینی ضریب ایمنی غارهای سنگی 
ارائه کردند. سپس یک آنالیز قابلیت اعتماد بر اساس معیار 
هوک براون برای پایداری جبهه کار تونل ارائه دادند. 

در این تحقیق سعی شده با بهره‌گیری از روابط 
هوک و براون و موهر کولمب» یک متدولوژی برای 
تحلیل پایداری دو تونل (در شرایط متفاوت) با استفاده 
از روش‌های قدرتمند قابلیت اعتماد شامل روش قابلیت 
اعتماد مرتبة اوّل. روش قابلیت اعتماد مرتبةٌ دوم و 
شبیه‌سازی مونت‌کارلو در نرم‌افزارهای جدید 1 و 
و۵8 ارائه شود. 


روش‌های ارزیابی قابلیت اعتماد 
با توجه به پیچیدگی محاسبهٌ انتگرال در رابطةٌ احتمال 
خرابی. محققان در پی یافتن روش‌هایی جهت تقریب 
زدن قابلیت اعتماد بوده‌اند. وظیفٌ روش‌های قابلیت 
عتماد محاسبةٌ احتمال هر رویداد مورد نظر است. هر کدام 
و وشن هابه شبزهای ابصمال شکیی وا ماه می که 
ز جمله روش‌های قابلیت اعتمادی که در این تحقیق 


ستفاده می‌شو ند می‌توان به روش مرتبهةٌ اول قابلیت 


عتماد» مرتبهةً دوم قابلیت اعتماد و روش شبه‌سازی 
مونت کارلو اشاره کرد. در ادامه به شرح این روش‌ها 


پرداخته می‌شود. 


نشریه مهندسی عمران فردوسی 


شاخص هاسفر -لیند (112:0107-100) و روش 

مرتبة اوّل قابلیت اعتماد 
شاخص هاسفر-لیند و يا شاخص قابلیت اعتماد. کمّیتی 
است که اجازه می‌دهد متغیرهای تصادفی نمود پیدا کنند 
و قضاوت در مورد پایداری بر اساس آنها انجام پذیرد. 
این اندیس یک مبنای مقایسه‌ای برای قابلیت اعتماد به 
شمار می‌رود و با استفاده از رابطةٌ زیر به دست می‌آید 
[14]: 


9 006-0 لجنه 8 


که در آن ۵6 بردار متغیرهای تصادفی (6) ۸ بردار 
میانگین متغفیرها. ‏ ماتریس کوواریانس و 3 دامن 
شکست است. در اینجا 8 از طریق یک مسألة بهینه‌سازی 
به دست خواهد آمد. لو و تنگ یک رابطٌ جایگزین برای 
رابطةٌ هاسوفر-لیند ارائه دادند که در آن به‌جای فضای 
استاندارد نرمال از فضای اصلی متغیرها استفاده می‌شود 
[14]: 


9 ۵ [۳ به ۱ 


که در آن ۶ ماتریس همبستگی و زه انحراف معیار 
متغیرهای تصادفی است. شاخحص قابلیت اعتماد را 
می‌توان به‌صورت هندسی نیز تعریف کرد. در فضای 
ترسیم یافته متغیرها به‌صورت استاندارد نرمال. شاخص 
قابلیت اعتماد را می‌توان فاصلهٌ بین مبداً و نقطه‌ای دانست 
کف نف ان نقطهٌ طراحی ()08ظ «عنوع) غیی گونته اک 
نقطه بر روی رویه حالت حدی قرار دارد و کمترین 
فاصله را تا هید دارد. یمن از به دست آمدان شاخضص 
قابلتتت؛ اعماد» استمال شکینت به.ظرنیق راظه. ابر 
محاسبه می‌شود [14]: 
63 (08 2۱ (9-8 < )۲ 

در اين رابطه ۵ بیانگر تابع توزیع تجمعی استاندارد 
ترمال است [14] 


روش مرتبةٌ دوم قابلیت اعتماد 
تابع حالت حدی ممکن است خطی يا غیر خطی باشد. 


نشريه مهندسی عمران فردوسی 


تحلیل پایداری یک تونل در شرایط هیدروستاتیک با اسفاده از... 


در روش مرتبةٌ اوّل. بسته به میزان غیر خطی بودن تابع 
فخاشتات: اال کیت داشار شعظا افو بصراکه 
تابع. با یک صفحه يا ابر صفحه تخمین زده می‌شود. 
زمانی که تابع حالت حدی انحنای زیادی داشته باشد. 
تقریب خطی, دفت محاسبات را پایین می‌آورد. در توابع 
غیر خطی روش مرتبةٌ دوم با کمک گرفتن از تقریب 
درجه دوم به یک جواب صحیح نزدیک می‌شود. این 
روش درواقع تکامل‌یافته روش مرتبه اوّل است و به‌جای 
استفاده از تقریب خطی, تابع حالت حدی را با یک 
!با له سس نی روف رک دیسا ۱۹۸۲ 
رابطهٌ زیر را برای محاسبهٌ احتمال شکست با تقریب مرتبة 
دوم ارائه داد [6]: ۱ 


۲ ۱۳ 
ره (ع۵ +۱ ۲1 (0-)۵ ه ۲ 
که در آن :1 نشان‌دهندة انحنای اصلی تابع حالت 


متغیر داشته باشد. (0-1) انحنا دارد [6], 


شبیه‌سازی مونت کارلو 
همان‌طور که از نام آن برمی‌آید کلیت‌های شبیه‌سازی 
مونت‌کارلو عبارت است از نمونه‌گیری به‌صورت تصادفی 
به‌منظور شبیه‌سازی تصنعی تعدادی زیادی آزمون و 
مشاهده نتایج آنها در حالت آنالیز قابلیت اعتماد سازه‌ها 
در ساده‌ترین شکل. احتمال خرابی در این روش به شکل 
زير محاسبه می‌شود [6]: 


(0 


در آنها ۰ > 6 است. تعداد آزمون‌ها آ به‌دقت مطلوب 
۴ ارتباط می‌یابد. روشن است که در روش مونت‌کارلو 
یک بازی شانس با توجه به خواص احتمالاتی معلوم 
به‌منظور حل مسأله به تعداد زیاد انجام می‌گیرد و سپس 
از این فرآیند نتیجه دلخواه یعنی احتمال خرابی استخراج 
می‌شود. از مزایای مهم اين روش عدم وجود هیچ 
محد ودیتی در نوع و شکل )1 و نیز تابع حالت حدی 
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ست [6]. 


بردارهای اهمَیّت 

بردارهای اهمیّت و حساسیت زیادی ارائه شده که 
اهمَیّت و تأثیرگذاری همة متفیرها را در قالب یک بردار 
بیان می‌کند و هر یک برای بیان این اهمّّت روش خاص 
خود را دارند که در اینجا به بردارهای اهمَیّت 0 و ۲۷ 


پرداخته می‌شود. 


بردار اهمیّت آلفا 
بردار اهمیّت آلفاء اهمَیّت هریک از متغیرها را بر اساس 
سهمی که در محاسبة واریانس تابع حالت حدی دارد بیان 
می‌کند. واریانس تابع حالت حدی ازاین‌جهت که در 
تایه تیال کته دا ایک بیان کته ابو 
اهمیّت است [16: 


۳۷ ۲ ۲ 1 ۲ 
۷]6[ < ۱۷۵۱ 0 ۷ ۹2 


سهم هر متغیر تصادفی در واریانس تابع حالت 
1 

حدی به‌اندازة 0 یعنی بردار یکه آن است. هرچه |زه| 
سل باشک اک ان ان که افعتت سار اس 
غاافت ۵ تد فر استمال شکست تاثر گذاز استر مقذاز 
مثبت » نشان می‌دهد با افزایش آن, متغیر ۲۶ افزایش پیدا 
می‌کند و مقدار منفی بیانگر آن است که با افزایش آن؛ 
متغیر ۳ کاهش می‌یابد [6]. 


برداز اهشتت انا 
این بردار همبستگی بین متغیرها را نیز در نظر می‌گیرد. 
به‌طوری‌که اگر همبستگی بین متغیرها صفر باشد بردار 
اهمیّت 0 و ۷ یکی هستند [6]. 
۸( 


گر ره ِ 
ارم ۷۲۰ 


که در آن 2 ماتریس انحراف معیار متغیرهای معادل 


نرمال و رو[ ماتریس ژاکوبین است. 


پایداری یک تونل در شرایط هیدروستاتیک 


وقتیکه در داخحل تودة سنگی که تحت تأثیرتنشیهای 
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برجا قرار دارد. تونل یا هر حفرةُ زیرزمینی دیگری 
احداث شود. وضعیت تنش موثر بر سنگ‌ها در مجاورت 
این تونل. به هم می‌خورد و توزیع جدیدی می‌یابد. برای 
تحلیل رفتار تونل» فرض می‌کنیم که تونل دایره‌ای شکلی 
به شعاع ,۲ تحت تنش هیدروستاتیک ,3 و نیز تحت‌فشار 
یکتوانعت دالعلی ۲ که سیستم نگهذاری توئل را اعمال 
می‌کند قرار دارد (شکل ۲ [6] 


(۹۱4۱۹۸۱ 
(۳۳9۳3۳9۳9۹9۹9 


زا 


شکل (۲): تتش موثر بر یک تونل دایره‌ای [15] 


اگر فشار داخلی مژثر بر تونل که ناشی از سیستم 
نگهداری نصب‌شده در آن است. از فشار بحرانی 
نگهداری ۳ کمتر شود در سنگ‌های اطراف تونل 
شکستگی رخ می‌دهد. فشار بحرانی م۲ از رابطة زیر به 


دست می‌آید [6]: 


0( < م۳ 


م6 مقاومت فشاری تک‌محوره سنگ است که با 


رابطة زیر تعریف می‌شود [6 


۲۷۰۵5 ۷ 


(۳( 


تحص 01 
۱-10 5 


نیز شیب منحنی تغییرات ,6 نسبت به 6 است 


و به‌صورت زیر بیان می‌شود (شکل ۳ [6]. 
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1 
3 
1 
9 
8 
دون < 
3 تنش محصور کننده 
شکل (۳): بیان ترسیمی معیار موهر کولمب [19] 
(م صنع+۱) 
(۱۱ چ << 1 
(م صنع) 


در اين روابط 6 تنش محوری که به ازای آن 
تکیت ی ۶ داده است. ,6 تنش محصورکننده. ) 
چسبندگی توده سنگ و ه زاویه اصطکاک داخلی سنگ 
است [6]. 

اک ققار دای ای کر او کار بحران و 
باشتله شگسگی ایساد ق شووی فا گنای اطر ات 
تونل الاستیک خواهد بود. در این صورت جابه‌جایی 
الاستیک شعاعی روبه درون دیواره تونل از رابطة زیر 
حاصل می‌شود [21) 


(«+۱) ۳ 
)۱۲ (ه- «) < یرل 


3 


که در آن ۷ ضریب پواسون و 5 مدول الاستیسیته 
است. در حالتی که فشار داخلی سیستم نگهداری ۲ کمتر 
از فشار بحرانی ۴ باشد. شعاع منطقة پلاستیک در 
اطراف تونل از رابطةٌ زیر به دست می‌آید [21]: 
۱ 
(- |سمج(اس ۱۳ 


)۱+۱()۲-۱ (۵+0 1" 
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در این صورت جابه‌جایی شعاعی کلی روبه درون 


دیوارة تونل خواهد شد [21]: 


۲ 


ِ ۱ بط - (۱-۷) ۲ حست < ,ولا 


1 ۳( " ب ۱( 


۱( 


معرفی نرم‌افزارهای به کارگرفته شده 
در این تحقیق از نرم‌افزارهای قابلیت اعتماد 1 و 
1 برای انجام تحلیل‌ها کمک گرفته شده است و در 
ادامه نتایج تحلیل‌ها توسط این دو نرم‌افزار در کنار هم 
آورده شده است. در این بخش به معرفی آنها پرداخته 


می‌شود. 


نرم‌افزار 11 
۲ نرم‌افزاری است که به‌وسیلاٌ آن می‌توان کارهای 
احتمالاتی و يا قابلیت اعتماد انجام داد. در این برنامه 
می‌توانیم پارامترهای مورد نظر را به‌صورت متغیر تصادفی 
به برنامه 1 داد. در این برنامه این امکان وجود دارد که 
برای تعداد زیادی متغیر ورودی این کار را انجام دهد و 
درنهایت میزان حساسیت هر یک از پارامترها را تعیین 
نماید. در این نرم‌افزار توزیع‌های مختلفی مانند نرمال» 
گامبل. بت نمایی» لاگ نرمال و رایلی وجود دارد که 
می‌توان برای برازش به یک سری داده از آنها استفاده 
نمود. انجام شبیه‌سازی مونت‌کارلو روش مرتبة اوّل 
قابلیت اعتماد و روش مرتبهة دوم قابلیت اعتماد از دیگر 
قابلیت‌های برنامه است. با انجام این آنالیزها می‌توان 
مشخص کرد برای یک تابع حالت حدی از یک سری 
متغیر تصادفی به ازای چه مقداری از آنها کمترین احتمال 
شکست وجود دارد. ضمناً این برنامه با برنامه‌هایی مانند 
طشه ۲5 همدص وعع‌قموم0۵ و ]۱۷۲۸ 
می‌تواند لینک شود. اين نرم‌افزار را می‌توان به‌صورت 
رایگان از طریق لینک (6/0ه ۵0ص اوننون) 


دریافت کرد. 
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نرم‌افزار »اهن018) 

پروژه‌ها به کار می‌رود که بر روی نرم‌افزار 0( 
17681 نصب می‌شود و به‌صورت یک نرم‌افزار افزودنی 
کار می‌کند. ۵1156 با تلفیق قابلیت‌های ۴۳6۵1 قادر به 
طراحی مدلی برای برآوردن نیازهای تحلیل ریسک با 
بهترین کاربری است. این نرم‌افزار برای تحلیل پروژه‌هایی 
با پارامترهای احتمالی بکار می‌رود. توزیم احتمال هر یک 
از پارامترهای مذکور بایستی قبلا مشخص شود. این 
پارامترها به همراه توزیعشان به‌عنوان ورودی وارد 
نرم‌افزار می‌شوند. سپس 6۵1856 با استفاده از روش‌های 
شبیه‌سازی مونت کارلو و با ۲3۷06۲0۷۲06 طناض1 و تولید 
اعداد تصادفی» مدل را بارها اجرا کرده و پارامترهای 
پروژه را تحلیل می‌کند و درنهایت مقادیر خروجی‌های 
مشخص‌شده و توزیم احتمال آنها را ارائه می‌کند. با 
پروژه را با توجه به انتظارات مورد نظر به دست آورد. 
این نرم‌افزار قادر است تا نتایج به‌دست‌آمده را به‌صورت 
گرافیکی نیز ارائه دهد که این قابلیت باعث فهم بهتر و 
سریع‌تر جواب‌ها می‌شود. تقریباً در تمام مدل‌هایی که با 
عدم قطعیت مواجه هستیم. می‌توان از اين نرم‌افزار 
استفاده کرد. 


تحلیل و نتایج 

در این بخش, نتایج به‌دست‌آمده از آنالیز قابلیت اعتماد 
مربوط به پایداری دو تونل ارائه شده است. تونل اول یک 
تونل دایره‌ای است که در یک فضای همگن همسان گرد 
و الاستیک حفرشده است و نحت‌فشار هیدروستاتیک ۳ 
و فشار نگه‌دارندة فاشلی یکنواخت ۳ قرار دارد [6]. 

تحلیل پایداری این تونل در سه نوع شرایط :۲ 
بررسی می‌شود که شامل مقادیر صفر ۰/۳ و ۰/۸ 
مگاپاسکال است و در هر بار بررسیء مقدار آن به‌صورت 
قطعی در رابطه‌ها وارد می‌شود. علاوه بر ز متغیرهای 
قطعی دیگر اين مسأله ۲.۵  <‏ 2۰.۳ ۸ ضریب 
پواسون و ۲ شعاع توئل هستند. همچنین متغیرهای 
تصادفی؛ 1۵ سین گی؛ 0 زاویة اصطکاکی داخلی و 1 
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تفای انا در جدول (۱) آورده شده است. هر سه پارامتر 
توزیع نرمال دارند. علاوه بر این متغیرهای ن) و 0 همیسته 
هستند و ضریب همبستگی آنها ۵- است. 


جدول (۱): مشخصات پارامترهای تصادفی تونل اول در توزیع 


نرمال 
متغیر تصادفی میانگین | انحراف معیار 
زاویهة اصطکاک داخلی (درجه) | ۲۲/۸۵ ۱/۳۱ 
چسبندگی (مگایاسکال) ۳۳ ۰/۸ 
مدول الاستیسیته (مگایاسکال) ۳۷۳ 1۸ 


پلاستیک و جابه‌جایی شعاعی کی به دست می‌آید. 
نخستین مرحله برای انجام تحلیل‌های قابلیت اعتماده 
تصمیم‌گیری در رابطه با معیارهای عملکرد به‌صورت 
ارائه یک تابع حالت حدی و پارامترهای مرتبط با آن 
است؛ بنابراین با استفاده از روابطی که گفته شد. توابع 
حالت حدی را به‌صورت زیر می‌سازيم. به‌طوری‌که اگر 
تونل از حد مجاز آنها تجاوز کند و بیشتر شود مقدار 
عددی به‌دست‌آمده برای تابع حالت حدی یک مقدار 
صورت مقدار عددی مثبت بوده و شرایط ایمن است. 


۳۳| 
(۱۵ 3 << 000 
)۱1( ت < ,ع < 0600 


در اين روابط با حد نهایی شعاع منطقه پلاستیک و 
۶ حد نهایی جابه‌جایی شعاعی تونل است. در اینجا 
مقدار ۳ < بآو مقدار ۰.۰۱ < ,5 است. به عبارتی تجاوز 
از این حد باعث خرابی تونل می‌شود و شکست در نقاطی 
اتفاق می‌افتد که مقدار تابع هالت شم مها راگن ستاونیک 
مقدار منفی شودد؛ بنابراین هدف, محاسبه احتمال شکست 
با استفاده از این توابع است. 


با استفاده از نرم‌افزار 137 اولین تحلیل را بر روی 
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تابع حالت حدی اوّل ((6,00) و ۰ < :3 انجام می‌دهیم؛ 
یعنی حالتی که تونل نگهداری نشده باشد. پس از اتمام 
آن به شکل () در صفحه ظاهر می‌شود. 

نرم‌افزار 11 برای رسیدن به مقدار شاخص قابلیت 
اعتماد از یک الگوریتم گام به گام استفاده می‌کند تا در 
فضای تبدیل يافته متغیرها به نقطه طراحی برسد. با توجه 
به حروجی که از نرم‌افزار 1*1 گرفته شده. نقطه‌ای که در 
فضای اصلی متغیرها دارای مختصات 
استاندارد نرمال دارای مختصات است. مولفهة اوّل در این 


اسنت: این نقطه در 
و در فضای 


مختصات مربوط به متغیر ن) و ملفه دوم مربوط به متغیر 
است. پس از مختصات نقطهٌ طراحی. بردار اهمَیّت 
گاما مشخص شده است. این بردار دارای دو ملفه است 
که ملفه اوّل به اهمَیّت متغیر ) و مژلفه دوم به اهمَیّت 
متغیر 0 می‌پردازد. در اینجا مقدار هر دو مژلفه منفی 
به‌دست‌آمده که این نشان‌دهندة آن است که با افزايش 
متغیرهای چسبندگی و زاویةٌ اصطکاک داخلی احتمال 
شکست کاهش می‌پابد. همچنین در اين بردان مقدار قدر 
مطلقی مولفة اوّل بیشتر از مولفةٌ دوم است که این نشان 


تحلیل پایداری یک توئل در شرایط هیدروستاتیک پا استفاده از... 


می‌دهد. متغیر چسبندگی نسبت به زاويةٌ اصطکاک داخلی 
بر فق, امال: شکب افر گلازی اش هرک 
تصمیمات مهندسی اولویّت بالاتری دارد. پس 
اهمیّت. آخرین نتایج. مربوط به شاخص قابلیت اعتماد و 
احتمال شکست است. در این تحلیل مقدار ۰.۷< 6 و 
۲ << :1 است. شکلی که نرم‌افزار 
تحروجی می‌دهد به‌صورت زير است. 
نموداری که در شکل (۵) مشخص شده نزدیک 
شدن گام‌های جستجو و کم شدن فاصلة آن به نقطة 
طراحی را نشان می‌دهد؛ بنابراین در گام آخر این فاصله 
نشان‌دهندة شاعص قابلیت. اعغماد است: در ادامه این 


از بردار 


۲ به‌عنوان 


تحلیل‌ها در مورد هر تونل. در هر دو نوع تابع حالت 
حدی» روش‌های قابلیت اعتماد به کارگرفته شده که نشان 
می‌دهد در مقدارهای مختلف از :۲ مقادیر متفاوتی از 
احتمال شکست به دست می‌آید. در این تحلیل‌ها از هر 
دو نرم‌افزار 1۳1و ۵1851 کمک گرفته شده و حروجی‌ها 
در قالب حداول (۲) و (۳) ارائه شده است. در حدول‌های 
زیر تحلیل‌هایی با تابع حالت حدی اوّل و دوم آورده شده 


است. 


)61 ۴۱۲1۵۱ ۷۷۲۲۳ ۵۵۸۷515 ۳0۴/۲ ۱6 ۲۵۲دک ربیب 


: 3 وع5۲ 


۳2610 < 2.38485 6-05 , ۳62 < 0۰000382833 , 091563۳766 ۶. 1 


۰ 20 2۲ 1 07۱۲ تیب 26 0۵1۳۲ 9۳ 51ع1 ع۲۳ ۳۵۲ 0۲ 52۳ 


053۳ 0۵1۳ 1۳ ۵۳191۳81 5866, 
[ ۵.185089 | 
۱ 0.401964 | 


۶ ,5۳862 53۳03۳0-۳0۳۵ ۷ کول وی 
۰50462 


-0 
| -۵.232502 | 


.۰ ۱2-35۸۳ پیت متس ۱ 
925ع- | 
۱ 32044 ۳ 


60۳۵۷۲1۳9 ۳ ۳5۵۳007 ۷۵۳ 1 291۶ ۷۰ 


۴۵۳۴۵۵۳ ۷۵۳12016 56۳5171۷117 25 2۳ 0۵۷۸۵ [۳ 


۰.۸۰۹ 


+ ۰ ۶ 


10 مخت ی دنرره 
31 2 86132 


:۳۵۶ ۲21 6۲ موی دزرره 
1 2 ۲۲ 


0۰ 391 0۰ 1۳ 0۵0۶ ۲515 لاد ۲06 جج یجید 


شکل (4): خروجی نرم‌افزار 1۲ برای تابع 0100 با روش قابلیت اعتماد مرتبةٌ اوّل 
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سال سی و سوم شماره دی ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی 


بح ع 2 9 سوت 


۷ 


و 


وان 


ح ۲ ۲ 
سپ م۱۷ وک 
شکل (۵): خروجی نرم‌افزار 181 برای تابع 0,0٩‏ با روش قابلیت اعتماد مرتبة اوّل 
جدول (۲): تخمین احتمال شکست در تابع حالت حدی اوّل 
نرم‌افزار 117 نرم‌افزار لهنع) 
فشار نگهداری (مگاپاسکال) روش مرتبة اوّل روش مرتبةٌُ دوم مونت‌کارلو مونت‌کارلو 
3 (۵: 3 (۳۶۵ (۳۶۵ (۳۶۵ 
۱/۷ ۳:۸۹ ۰/6 ۳9/۸۷۷ ۳۳/۸/۹ ۳۶ 
۳ ۳/۲ 1/۰۹ ۶ 1/۰۹ 1/۰۹ 
۰/۸ ۵ | ۶/۷۱۰۲ ۷/1۵ 2 
جدول (۲): تخمین احتمال شکست در تابع حالت حدی دوم 
نرم‌افزار 117 نرم‌افزار لهن8) 
فشار نگهداری (مگاپاسکال) روش مرتبة اوّل روش مرتبةٌ دوم مونت‌کارلو مونت‌کارلو 
3 (۳۶۵ 3 (۳۶۵ (۵ (۳۶۵ 
۰/۸ ۱/۷ ۰/۹ ۱/۷ ۰/۹۵ ۱,۹۲ ۱/۸ 
در اینجا تعداد اعداد تصادفی که در شبیه‌سازی احتمال شکست دارد. 


مونت کارلو تولید می‌شود ‏ ۱۰۰۰ < 3۷ است. شکل (1) 
روند همگرایی روش شبیه‌سازی مونت‌کارلو به جواب را 
نشان می‌دهد. با توجه به این شکل. تغییرات احتمال 
شکست» حتی تا آخرین دادهُ شبیه‌سازی‌شده در حال 
کاهش است و منجر به یک عدد ثابت نشده است. 
می‌توان گفت این تغییرات در محاسبه احتمال شکست. 
سکن ابیت: با ادامه شببه‌سازی, کاهین. بدا کند.. النتد 


تصمیم به ادامة شبیه‌سازی بستگی به‌دقّت قابل قبول 


سال سی و سوم شماره دی ۱۳۹۹ 


نرم‌افزار 61851 برای محاسبةٌ احتمال خرابی از 
روش مونت‌کارلو استفاده می‌کند. در شکل زیر خروجی 
ان نرم‌افزار برای تابع ان عفن این | مسوف نت ات 
که در آن بر اساس مساحت محصورشده در زیر تابع 
چگالی احتمال ()) مقدار احتمال خرابی به دست 
می‌آید. احتمالی که در شکل (۷) به‌دست‌آمده ۲۶/ است 
که تقریباً به جوابی که با نرم‌افزار 1 به دست آمد نزدیک 


انتتتتا: 


نشریه مهندسی عمران فردوسی 


1۸ 


با توجه به اعداد نمایش داده شده در جداول (۲) و 
(۳) احتمال شکست‌ها در هر روش و برای هر تابع حالت 
حدی متفاوت است که از این احتلاف در احتمالات 
برداشت‌هایی می‌توان کرد. زمانی که مقدار فشار 
نگهداری صفر باشد و به عبارتی شرایط بدون نگهداری 
فرض شود. اتتمال شکست. مجاسیه 1 : ( ۲۶ تیاو 
بالا خواهد بود؛ بنابراین تونل مورد نظر با وضعیتی که در 
مورد سایر پارامترها اعلام شده ایمن نیست و نیاز به یک 


2 تا << 


7297.2973, 8 


تحلیل پایداری یک توئل در شرایط هیدروستاتیک با استفاده از... 


سیستم نگهداری با فشار مناسب دارد. از طرفی احتمال 
شکست. زمانی که یک سیستم نگهداری با فشار 
محافظت کننده ۰/۳ اعمال شود بسیار کوچک (۰/۰5) 
برآورد شده است؛ بنابراین می‌تواند برای ایمن‌سازی تونل 
مورد نظر مناسب باشد. اگرچه فشار نگهداری ۰/۸ نیز 
شرایط ایمنی را فراهم می‌کند اما با توجه به هزينة بالایی 
که به همراه دارد مناسب نخواهد بود. 


صهدهاه ۲ عطا ۵۴ ام 
3 


۰ 
خک 
4 


لت 
رح 


شکل (1): روند همگرایی شبیه‌سازی مونت‌کارلو به احتمال شکست نهایی 


4 ی‌تاوتاها5 
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47/553 رت 
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مججرو ۱۵ 
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بو ۱۵0 
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1۰ 22/433 ۹ 
۱ ۱۷۵0 
موچو ۱۷۵4 
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وو- ۷ ۱ 
0/04 2 ۱ ۳ 3 ضُ ه‌ 
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شکل (۷): خروجی نرم‌افزار ۵18 برای تابع 0100 
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هادی فتاحی- فاطمه جیریایی شراهی 


در مقادیر بالای احتمال شکست. یعنی در تونل اوّل 
با فشار نگهداری صفر. روش قابلیت اعتماد مرتبهٌ اوّل و 
روش شبیه‌سازی مونت‌کارلو خروجی‌های نزدیک به هم 
دارند. اما این خروجی‌ها نسبت به احتمال شکست روش 
قابلیت اعتماد مرتبةً دوم کمی اختلاف دارد. این اختلاف 
در مقادیر احتمالات ممکن است به دلیل نزدیکی رویه 
حالت حدی اول ۰ < () ,0 به یک خط باشد که 
روش قابلیت اعتماد مرتبةٌ دوم را دجار حطا می کند. 
چراکه در روش مرتبةٌ دوم از یک سهمی برای تقریب‌زدن 
رویهةٌ حالت حدی استفاده می‌شود. از طرفی در احتمال 
شکست‌های خیلی کوحک. سه روش یادشده خروجی 
شبیه به هم دارند و از آنجا که شبیه‌سازی مونت‌کارلو 
برای رسیدن به دقتی بالا در احتمالات بسیار کوچک. نیاز 
به تعداد زیادی شبیه‌سازی و صرف زمان زیادی دارد. 
استفاده از آن در این شرایط توصیّه نمی‌شود؛ بنابراین با 
توجه به پیچیدگی تابع حالت حدی دوم و دقت بالای 
برای روش شبیه‌سازی مونت‌کارلو باشد. 

در هر دو تونل برای تخمین احتمال شکست‌ها از 
دو تابع حالت حدی ( 6 و 00 استفاده شد که بر 
اساس نتایجی که به دست آمد در هر حالت خروجی این 
دو تابع با هم متفاوت بوده و در همه حالت‌ها احتمال 
شکست به‌دست‌آمده از (6,)6 نسبت به (6,)06 بیشتر 
بوده است. این ممکن است به دلیل در نظر گرفتن متغیر 
در 6,600 باشد که از ان نظر می‌تواند تخمین واقمیتر 
و دقیق‌تری از احتمالات داشته باشد. 

تونل دوم نیز همانند تونل اوّل یک تونل دایره‌ای 
است که در یک فضای همگن. همسانگرد و الاستیک 
حفرشده آنحتت و تحت‌فشار هیدروستاتیک ۰ و فشار 
نگه‌دارنده داخلی یکنواخت ‏ قرار دارد [19]. متغیرهای 
قطعی این مسأله ۲.۵ < .۳ ۰.۳< ۷ ,۲ شعاع تونل 
هستند. همچنین متغیرهای تصادفی» 0 یبتک کیان 0 
زاویه اصطکاک داخلی» تل مدول الاستیسیته و ۳ فشار 


نگه‌دارنده داخلی هستند که نوزیع» میانگین و انحراف 


سال سی و سوم شماره دی ۱۳۹۹ 


1۹ 


معیار آنها در جدول (4) آورده شده است. هر سه پارامتر 
توزیع نرمال دارند. علاوه بر این متغیرهای ) و 0 همبسته 
هستند و ضریب همبستگی آنها 0 اسستا: 


جدول (4): مشخصات پارامترهای تصادفی تونل دوم 


متغیر تصادفی | تابع توزیع | میانگین | انحراف معیار 
زاویة اصطکای 
نرمال ‏ | ۲۱/۰۹۳ ۱۳/۳1۵ 
داخلی (درجه) 
5 
چسبن کی نرمال ۰۳۸۷ ۰۳۸ 
(مکاپاسکال) 
فشار نگهداری 
تیان / ۰/۹ 
(مکاپاسکال) 
مدول الاستیسیته 
نرمال ۳۷۳ 1۸ 
(مگایاسکال) 


همچنین در این مسأله مقدار 2حاو 2۰,۱۲ ره 
ات تحلیل‌هایی که انجام شده به‌طور خللاصه در جدول 


(6 آورده شتلاه: اسست: 


جدول (۵): تخمین احتمال شکست در تونل دوم 


روش مرتبةٌ دوم 
(1) 


4 


600 | 62060 | 60۵ | )62 
احتمال شکست (1) 1/۵ ۳/۵ 1/۹ 3 


در مورد اين مسأله فشار نگهداری به‌عنوان یک 
متغیر تصادفی با یک میانگین و انحراف معیار ارائه شده 
است. این متغیر با در نظر گرفتن باه اطمینان (۳5 ع 6 
می‌تواند بین مقادیر (۰/۲۲-۰/۵۸) تغییر کند. با مقايسةٌ 
این مقادیر نسبت به مقادیر در نظر گرفته شده برای فشار 
نگهداری در تونل اوّل. می‌توان دریافت تونل اوّل به سبب 
تفاوت در مقادیر پارامترهای تصادفی آن. نسبت به تونل 
دوم ایمن‌تر بوده و فشار نگهداری مورد نیاز برای رسیدن 
به احتمال شکست قابل قبول برای آن کمتر بود است. 


شکل‌های (۸) و )٩(‏ نمودارهایی را نمایش می‌دهد که 
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بمکو نکن تیان شاتفشن نی اعشان .ی اسیال 
تکیت وا قیال هی تفای بکدا ره اسان 
سفق مان طرن. که ین سنا آفایی قفا 
نگهداری, نمودار تغییرات احتمال شکست با شیب تندی 
کاهش می‌یابد و تقریباً در مقدار ۰.۲ < مگاپاسکال 
این نمودار ثابت می‌شود. معمولاً احتمال شکست قابل 


قبول فرافه وهای ای ۵ ام کهآ گر کین 
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0۳25 
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تحلیل پایداری یک تونل در شرایط هیدروستاتیک با اسفاده از... 


باشد با تدارک یک سیستم نگهداری با فشار نگهداری 
تقریا ۵/۱۲ مکاباشکال: اما فنکست به معدان فورد 
نظر کاهش می‌یابد. همچنین نمودار شاحص قابلیت 
اعتماد با افزایش 8 با یک شیب ثابت افزایش می‌یابده 
به‌طوری که افزایش شاخصی قابلیت اعتماد کاهش احمال 
شکست را نشان می‌دهد. 


1 08 06 


فشار نگهداری مگاپاسکال 


شکل (۸): تأثیر فشار نگهداری بر احتمال شکست 


۳ 
۵ دك ره 0 0 ۱۵ 
اعتماد 


‌ِ 
0 
ت‌ِ 


1 08 06 


فشار نگهداری مگاپاسکال 


شکا (64: تاثیر فشار نگهٌذاوی بر شناغتصن قانلیت اغتماد 
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هادی فتاحی- فاطمه جیریایی شراهی 


با استفاده از تحلیل حساسیتی که بر مبنای آنالیز 
قابلیت اعتماد انجام شده است. نتایجی به این شرح به 
دست آمده است. آنالیز تحلیل حساسیت که به‌منظور 
تعیین اهمَیّت هر یک از پارامترها و همچنین نقش آنها 
در سیستم و احتمال شکست آن می‌پردازد. در 
نرم‌افزارهای مختلف ممکن است توسط روش‌های 
متفاوتی انجام شود. در نرم‌افزار 1*۲ به دلیل وجود 
همبستگی بین متغیرها از بردار اهمَیّت گاما استفاده شده 
است. در شکل‌های (۱۰) و (۱۱) نمودارهایی جهت 
نمایش اهمَیّت متغیرها در هر دو تونل اوّل و دوم و اعداد 
آن به‌طور واضح‌تر در جداول (1) و (۷) نشان داده شده 
است. به‌طوری که اهمیّت متغیرها با توجه به نوع برنامة 
استفاده شده و تابع حالت حدی مربوطه مشخص شده 
است. 

با توجه به اين که تابع حالت حدی دوم دارای متغیر 


۲ نیز است. نمایش بهتری را از همه متغیرهای تصادفی 


12 
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10۱5۱/۶102( 


11/0100 


15/۶202( 
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ارائه می‌دهد. بر اساس این نتایج در تونل اوّل وقتی از 
(06) ,6 استفاده می‌شود متغیر چسبندگی دارای اهمیّت 
بالاتری است و زمانی که از (6,)28 استفاده می‌شود با 
اضافه شدن متغیر 8 از دو متغیر قبلی پیشی می‌گیرد و 
دارای اهمَیّت بیشتری می‌شود. همچنین در مورد تونل 
دوم مدول الاستیسیته دارای اهمّیّت بالایی است؛ بنابراین 
می‌توان گفت مدول الاستیسیته نسبت به سایر متغیرها در 
طراحی بهینه سیستم نقش مهمی دارد. در مورد سایر 
مرها آهشکه آنها دز هر یی از مالت‌ها مضاویخ اسست: 
برای مثال در مورد تونل اول بعد از _ ) یعنی چسبندگی 
اهمَیّت بیشتری دارد؛ اما در مورد تونل دوم زاوية 
اصطکاک داخلی جایگاه دوم در مهم‌ترین متغیرها را دارد؛ 
بنابراین ممکن است برای هر مطالعةٌ موردی یک تحلیل 


حساسیت ویژه مطرح شود. 


۶" 
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شکل (۱۰ تحلیل حساسیث بر مبنای آنالیر قابلیت اعتماد در توتل ال 
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تحلیل پایداری یک توئل در شرایط هیدروستاتیک با استفاده از... 
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شکل (۱۱): تحلیل حساسیت بر مبنای آنالیز قابلیت اعتماد در تونل دوم 


جدول (1): تحلیل حساسیت بر مبنای آنالیز قابلیت اعتماد در تونل 


نرم‌افزار/تابع حالت حدی چسبندگی 
6۵8/0۱00 ۹۷ 
81/0۱00 ۲ 
689/0200 ۰/۳۶ 
81/0200 ۳۲ 


ول 
میزان اهمَیّت 
زاوية اصطکاک داخلی لول الاسکیسه 
۳۹ ۰ 
۰۳۸ 
۳۲ ۴۵ 
۳۳ ۹۲ 


جدول (۷): تحلیل حساسیت بر مبنای آنالیز قابلیت اعتماد در تونل دوم 


میزان اهمَیّت 
نرم‌افزار/تابع حالت حدی چسبند گی زاوية اصطکاک داخلی | فشار نگهداری (۷]۳2) | مدول الاستیسیته 
6۵۳9/0100 ۰/۳۰ ۰:0۵ 1۲ 
۳/0100 ۹ ۰/۷۷ ۰/۰ 2 
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نتیجه گیری نشان داد در مورد تابع حالت حدی اوّل در احتمال 


مرتبةٌ دوم و شبیه‌سازی مونت‌کارلو برای تحلیل پایداری 


نشريه مهندسی عمران فردوسی 


صفر). روش قابلیت اعتماد مرتبةٌ اوّل و شبیه‌سازی 
مونت‌کارلو تقریباً شبیه به هم عمل می‌کنند. از طرفی در 


سال سی و سوم شماره دی ۱۳۹۹ 


هادی فتاحی- فاطمه جیریایی شراهی 


احتمال شکست‌های خیلی کوچک» هر سه روش 
خروجی نزدیک به هم دارند. البته شبیه‌سازی مونت کارلو 
در مقادیر احتمال بسیار کوچک و دقت بسیار بالاء نیازمند 
تولید تعداد زیادی داد رندم و شبیه‌سازی هریک از آنها 
است و هزین محاسباتی بالایی دارد. از اين نظر استفاده 
از آن در این شرایط توصیّه نمی‌شود؛ بنابراین روش 
قابلیت اعتماد مرتبةٌ دوم با توجه به‌دقت بالای آن در 
تقریب زدن توابع» با وجود پیچیدگی تابع حالت حدی 
دوم. می‌تواند جایگزین مناسبی برای روش شبیه‌سازی 


۳ 


همچنین نتایج تحلیل حساسیت بر روی متغیرها 
تشان داد که مدول الا مشستیه دارآ هت وتات کار 
بسیار بالایی است و در تحلیل‌ها به لحاظ اهمَیّت و 
تأثیرگذاری اولویّت اوّل را دارد. به‌علاوه نتایج نشان داد 
که شاخص قابلیت اعتماد وابستگی زیادی به فشار 
نگه‌دارنده دارد. تونل‌های مورد بررسیء در حالت بدون 
برای آنها سیستم نگهداری تعریف شود که اگر در تونل» 
فشار سیستم نگهداری ۰.۳ تِ ۳ مکاپاسکال ایجاد شود. 


کل مت اختمال شکسست: قایز: قبولین خخواهد داششت: 
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11 تحلیل پایداری یک توثل در شرایط هیدروستائیک با اسفاده از... 
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